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ABSTRACT 

Using previous measurements and quantum chemical calculations to derive 
the molecular properties of the TiH molecule, we obtain new values for its ro- 
vibrational constants, thermochemical data, spectral line lists, line strengths, and 
absorption opacities. Furthermore, we calculate the abundance of TiH in M and L 
dwarf atmospheres and conclude that it is much higher than previously thought. 
We find that the TiH/TiO ratio increases strongly with decreasing metallicity, 
and at high temperatures can exceed unity. We suggest that, particularly for 
subdwarf L and M dwarfs, spectral features of TiH near ~0.52 /iin, 0.94 /im, 
and in the H band may be more easily measurable than heretofore thought. 
The recent possible identification in the L subdwarf 2MASS J0532 of the 0.94 
jttm feature of TiH is in keeping with this expectation. We speculate that looking 
for TiH in other dwarfs and subdwarfs will shed light on the distinctive titanium 
chemistry of the atmospheres of substellar-mass objects and the dimmest stars. 

Subject headings: infrared: stars — stars: fundamental parameters — stars: low 
mass, brown dwarfs, subdwarfs, spectroscopy, atmospheres, spectral synthesis 
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1. INTRODUCTION 

The discovery of the new spectroscopic classes L and T (Nakajima et al. 1995; Oppen- 
heimer et al. 1995; Burgasser et al. 2000abc; Kirkpatrick et al. 1999) at lower effective 
temperatures, (T eff ), (2200 K^700 K) than those encountered in M dwarfs has necessi- 
tated the calculation (and recalculation) of thermochemical and spectroscopic data for the 
molecules that predominate in such cool atmospheres. At temperatures (T) below 2500 K, 
a variety of molecules not dominant in traditional stellar atmospheres become important. 
Though a few molecules, such as TiO, VO, CO, and H 2 0, have been features of M dwarf 
studies for some time, for the study of L and T dwarfs, even for solar metallicities, the 
molecules CH 4 , FeH, CrH, CaH, and MgH take on increasingly important roles. The last 
four metal hydrides are particularly important in subdwarfs, for which the metallicity is 
significantly sub-solar. In subdwarf M dwarfs, the spectral signatures of "bimetallic" species 
such as TiO and VO weaken, while those of the "monometallic" hydrides, such as FeH, 
CaH, CrH, and MgH, increase in relative strength. This shift from metal oxides to metal 
hydrides in subdwarfs is an old story for M dwarfs, but recently, Burgasser et al. (2003) 
and Burgasser, Kirkpatrick & Lepine (2004) have discovered two L dwarf subdwarfs, 2MASS 
J05325346+8246465 and 2MASS J16262034+3925190, and this shift to hydrides is clear in 
them as well. In fact, Burgasser, Kirkpatrick & Lepine (2004) have recently tentatively 
identified titanium hydride (TiH) in 2MASS J0532 at ~0.94 /xm, a molecule that heretofore 
had not been seen in "stellar" atmospheres (see also Burgasser 2004). 

In support of L and T dwarf studies, we have an ongoing project to calculate line lists, 
line strengths, abundances, and opacities of metal hydrides. Burrows et al. (2002) and 
Dulick et al. (2003) have already published new calculations for CrH and FeH, respectively. 
With this paper, and in response to the recent detection of TiH in 2MASS J0532, we add 
to this list new spectroscopic, thermochemical, and absorption opacity calculations for TiH. 
In §2, we summarize the spectroscopic measurements we use to help constrain the spectro- 
scopic constants of TiH. We continue in §3 with a discussion of the computational chemistry 
calculations we performed in conjunction with our analysis of the laboratory work on TiH. 
In §5, we calculate TiH abundances with the newly-derived thermochemical data and par- 
tition functions and find that the TiH abundances, while not high, are much higher than 
previous estimates. We also describe the updates to our general chemical abundance code 
that are most relevant to titanium chemistry. Section 6 describes how we use these lists to 
derive opacities at any temperature and pressure and §7 summarizes our general conclusions 
concerning these updated TiH ro-vibrational constants, abundances, and opacities. 
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2. SUMMARY OF SPECTROSCOPIC WORK 

Although spectra of TiH have been available for several decades (Smith & Gaydon 1971), 
major progress on the spectral assignment has been made only recently (Steimle et al. 1991; 
Launila & Lindgren 1996; Andersson, Balfour & Lindgren 2003; Andersson et al. 2003). The 
spectrum of TiH was first observed by Smith and Gaydon in 1971 in a shock tube and a 
tentative 4 A — 4 $ assignment was proposed for a band near 530 nm (Smith & Gaydon 1971). 
By comparing the TiH measurements (Smith & Gaydon 1971) with the spectra of a Ori, a 
Sco, and 5 Vir, Yerle (1979) proposed that TiH was present in stellar photospheres. Yerle's 
tentative TiH assignments in these complex stellar spectra are very dubious. 

The first high-resolution investigation of TiH was made by Steimle et al. (1991) who 
observed the laser excitation spectrum of the 530-nm band. After rotational analysis, this 
band was assigned as the 4 r 5 / 2 — Y 4 $3/ 2 sub-band of a 4 r — X 4 $ transition (called the 
B 4 T — X 4 $ transition in the present paper). The spectrum of the 530-nm transition of 
TiH was later investigated by Launila and Lindgren (1996) by heating titanium metal in an 
atmosphere of about 250 Torr of hydrogen, and the spectra were recorded using a Fourier 
transform spectrometer. The assignment of the 530-nm system was confirmed as a 4 r — X 4 $ 
transition by the rotational analysis of the four sub-bands of the 0-0 vibrational band. The 
spectrum of the analogous 4 r — X 4 $ transition of TiD has recently been studied at high 
resolution using a Fourier transform spectrometer (Andersson, Balfour & Lindgren 2003). 
TiD was also produced by heating titanium powder and 250 Torr of D 2 in a King furnace. 
The rotational analysis of the ^7/2 — Y 4 $ 5 / 2 and 4 r 5 / 2 — Y 4 $ 3 / 2 sub-bands produced the 
first rotational constants for TiD (Andersson, Balfour & Lindgren 2003). 

In a more recent study, the spectra of TiH and TiD were observed in the near infrared 
(Andersson et al. 2003) using the Fourier transform spectrometer of the National Solar Ob- 
servatory at the Kitt Peak. The molecules were produced in a titanium hollow cathode 
lamp by discharging a mixture of Ne and H 2 or D 2 gases. A new transition with complex 
rotational structure was observed near 938 nm and was assigned as a 4 $ — X 4 $ transition 
(called the A 4 $ — X 4( 3> transition in the present paper). The complexity of this transition is 
due to the presence of perturbations, as well as overlapping from another unassigned transi- 
tion (Andersson et al. 2003). The spectroscopic constants for the TiH and TiD states were 
obtained by rotational analysis of the 0-0 band of the four sub-bands. 

In other experimental studies of TiH, a dissociation energy of 48.9±2.1 kcal mole -1 was 
obtained by Chen et al. (1991) using guided ion-beam mass spectrometry. The fundamental 
vibrational intervals of 1385.3 cm -1 for TiH and 1003.6 cm -1 for TiD were measured from 
matrix infrared absorption spectra (Chertihin & Andrews 1994). The molecules were formed 
by the reaction of laser-ablated Ti atoms with H 2 or D 2 and were isolated in an argon matrix 
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at 10 K. 

There have been several previous theoretical investigations of TiH. Ab initio calculations 
for TiH and other first-row transition metal hydrides have been carried out by Walch and 
Bauschlicher (1983), Chong et al. (1986), and Anglada et al. (1990). Spectroscopic prop- 
erties of the ground states of the first row transition metal hydrides were predicted in these 
studies. Chong et al. (1986) and Anglada et al. (1990) have also calculated spectroscopic 
properties of some low-lying states. In a recent study the spectroscopic properties of the 
ground and some low-lying electronic states of TiH were computed by Koseki et al. (2002) 
by high level ab initio calculations that included spin-orbit coupling. The potential energy 
curves of low-lying states were calculated using both effective core potentials and all-electron 
approaches. 

3. COMPUTATIONAL CHEMISTRY METHODS AND RESULTS 

To supplement the sparse spectroscopic measurements summarized in §2, we calculate 
the spectroscopic constants, transition dipole moments, line strengths, and ro- vibrational 
constants of the TiH molecule using standard ab initio quantum chemical techniques. The 
orbitals of the TiH molecule are optimized using a state-averaged complete-active-space self- 
consistent-field (SA-CASSCF) approach (Roos, Taylor & Siegbahn 1980), with symmetry 
and equivalence restrictions imposed. Our initial choice for the active space include the Ti 
3d, 4s, and 4p orbitals and the H Is orbital. In Civ symmetry, this active space corresponds to 
four ai, two &i, two 62, and one 02 orbital, which is denoted (4221). Test calculations showed 
that this active space is much too small to study both the B A T— X 4 $ and A 4 $— X 4 $ systems 
simultaneously. In fact, no practical choice for the active space was found that allowed the 
study of both band systems simultaneously. Therefore, the two band systems are studied 
using separate calculations. For the _B 4 r — X 4 <3> system, the (4221) active space is sufficient. 
In this series of calculations, two 4 A states and one state each of 4 n, 4 $, and 4 r are included 
in the SA procedure, with equal weight given to each state. In the A 4 $ — X 4 $ SA-CASSCF 
calculations, it was necessary to increase the active space, and our final choice is a (7222) 
active space that has two additional a orbitals and one additional S orbital, relative to the 
(4221) active space. Only the two 4 $ states are included in the SA-CASSCF calculations 
for the A 4 $ - X 4 $ system. 

More extensive correlation is included using the internally contracted multireference con- 
figuration interaction (IC-MRCI) approach (Werner & Knowles 1988; Knowles & Werner 1988). 
All configurations from the CASSCF calculations are used as references in the IC-MRCI cal- 
culations. The valence calculations correlate the Ti 3d and 4s electrons and the H electron, 
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while in the core+valence calculations the Ti 3s and 3p electrons are also correlated. The Ti 
3s- and 3p-like orbitals are in the inactive space so they are doubly occupied in all reference 
configurations. The effect of higher excitations is accounted for using the multireference 
analogue of the Davidson correction (+Q) (Langhoff & Davidson 1974). Since the IC-MRCI 
calculations are performed in Ci v symmetry, the A and T states are in the same symmetry as 
are the $ and II states. Thus, for the B 4 T — X 4 $ calculation, the IC-MRCI reference space 
included three roots of 4 Ai symmetry and two roots of 4 Bi symmetry. For the A 4 $ — X 4 $ 
system, the reference space includes four roots of 4 Bi symmetry since there are two 4 II states 
below the A 4 $ state. 

Scalar relativistic effects are accounted for using the Douglas-Kroll-Hess (DKH) ap- 
proach (Hess 1986). For Ti we use the (21sl6p9d4f3glh)/[7s8p6d4f3glh] quadruple zeta 
(QZ) 3s3p basis set (Bauschlicher 1999), while for H we use the correlation-consistent po- 
larized valence QZ set (Dunning 1989). The nonrelativistic contraction coefficients are 
replaced with those from DKH calculations. All calculations are performed with Mol- 
pro2002 (Werner & Knowles 2002) that is modified to compute the DKH integrals [C.W. 
Bauschlicher, unpublished] . 

The transition dipole moments for the forbidden 4 n- 4 r and 4 n- 4 <3> transitions are close to 
zero. Therefore, we conclude that the states are quite pure despite the fact that the IC-MRCI 
calculations are performed in Civ symmetry. We summarize the computed spectroscopic con- 
stants along with available experimental data (Andersson et al. 2003; Steimle et al. 1991; 
Launila & Lindgren 1996; Chertihin & Andrews 1994) in Table 1. We do not include the 
(2) 4 A state from the B 4 T -X 4 $ calculation, since, based on the weight of the reference con- 
figurations in IC-MRCI calculations, this state is not well described by these calculations. 

We first compare the results for the X 4 $ state as a function of treatment. For the 
treatments that include only valence correlation, A 4 $ — X 4 Q and _B 4 r — X 4 $ calculations 
yield very similar results for the X 4 $ state. When core correlation is included, there is 
a larger difference between the two sets of calculations. For example, note that the u e 
values (see Table 1) which differ by 6.4 cm -1 at the valence level, differ by 19.7 cm -1 
when core correlation is included. Our computed u e values are significantly larger than 
the value assigned by Chertihin and Andrews (1994) to TiH in an Ar matrix. However, 
recently L. Andrews [personal communication] has suggested that their value is incorrect. 
The A 4 $ — X 4 $ core+valence treatment results in a small increase in the dipole moment, 
which is now in excellent agreement with experiment, while for the -B 4 r — X 4 $ treatment, 
the inclusion of core correlation leads to a significant decrease in size and a value that 
is significantly smaller than experiment. Using a finite field approach for the X state in 
the _B 4 r — X 4 $ core+valence treatment yields a value that is in better agreement with 
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experiment, but is still 0.29 D too small. The A 4 & — X 4 & core+valence treatment yields the 
best agreement with experiment for the value of the Ti-H bond distance (r ) of the X state, 
in addition to yielding the best dipole moment. 

The core+valence treatment yields a better r value for the A 4 & state, but the dis- 
agreement with experiment is larger than for the X state. The valence treatment yields a 
somewhat better term value (T ), with both the valence and core+valence values being larger 
than experiment. The inclusion of core correlation improves the agreement with experiment 
for the r value for the B 4 T state and also improves the agreement with the experimental 
To. While the change in most uo e values with the inclusion of core correlation is small, the 
inclusion of core correlation leads to a sizable reduction in u e for the A 4 Q state. We observe 
a decrease in the dipole moments of the 4 n, 4 A, and B 4 T states when core correlation is 
included, as found for the X state in the B 4 T — X 4 $ treatment. As for the X state, the va- 
lence level dipole moment for the B 4 T state agrees better with experiment. We suspect that 
true values for the 4 n and 4 A states are closer to those obtained in the valence treatment 
than those obtained in the core+valence treatment. 

While the dipole moments vary significantly when core correlation is included, the effect 
of including core correlation on the dipole transition moment is only up to about 3 percent 
in the Franck-Condon region for both transitions of interest. For the A 4 $ — X 4 $ transition, 
core correlation increases the moment, while for the _B 4 r — X 4 $ transition core correlation 
decreases the moment. 

Given the limited experimental data it is difficult to definitively determine the best 
values for u e , u e x e , and the transition moment. On the basis of previous calculation on CrH 
(Bauschlicher et al. 2001; Burrows et al. 2002) and FeH (Dulick et al. 2003), we feel that the 
calculations including the core+valence may yield superior values for all spectroscopic con- 
stants excluding the dipole moment, where the smaller active space used in the -B 4 r — X 4 $ 
calculations results in the valence treatment being superior for this property. For the X 4 $ 
state, we believe the results obtained using the A 4 & — X 4 $ treatment are superior because 
it has a larger active space and fewer states in the SA procedure. We also note that for 
the A 4 $ — X 4 $ treatment, the X state dipole moment improves when core correlation is 
included, which also supports using the A 4 § — X 4 $ treatment for the X state. In the end, 
transition dipole moment functions computed at the core+valence level were used to gen- 
erate the Einstein A v / V » values for the -B 4 r — X 4 $ and A 4 $ — X 4 $ transitions in Table 2. 
The transition moments and IC-MRCI+Q potentials used to compute the A values in Table 
2 are plotted in Figures 1 and 2, respectively. 

Finally, we compare our best values to previous theoretical values. Our best results are 
in reasonable agreement with those reported by Anglada et al. (1990); for the X and A 
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states their bond distances are 0.07 and 0.02 A longer than our best values, respectively, 
and therefore in slightly worse agreement with experiment. Their X state u e value is about 
50 cm -1 smaller than our best value, while their A state value is about 140 cm' 1 larger than 
our best value. For the one bond distance for which they report their transition moments, 
their A — X transition moment differs from our value by about 10%, while their B — X 
moment is only very slightly smaller than our value. The X state r e values reported by 
Koseki et al. (2002) are about 0.1 A longer than our r e value. Their all electron oo e value is 
in better agreement with our value than is their value with scalar relativistic effects added. 
This is a bit surprising since our calculations include the scalar relativistic effects. Overall 
the agreement between the different theoretical approaches is relatively good. Our use of 
much larger basis sets and the inclusion of core-valence correlation results in bond lengths 
that are in better agreement with experiment, and we assume that the other spectroscopic 
properties are more accurate in our calculations as well. 



4. LINE POSITIONS AND INTENSITIES 

A common feature associated with the electronic structure of transition metal hydrides 
is that the excited electronic states are often heavily perturbed by unobserved neighboring 
states whereas the ground states are for the most part unperturbed. The X 4 $, A 4( 3>, and 
B 4 T states of TiH are no exception. Recent work on the analyses of high-resolution Fourier 
transform emission and laser excitation spectra of the A 4 $ — X 4 <£> (Andersson et al. 2003) 
and B 4 T - X 4 $ (Steimle et al. 1991; Launila & Lindgren 1996) - bands have finally 
yielded rotational assignments for the heavily perturbed A and B states. 

As a prerequisite to calculating opacities for temperatures up to 2000 K, a synthesized 
spectrum was generated comprised of rovibrational line positions and Einstein A-values for 
quantum numbers up to v" = v' = 5 and J" = J' = 50.5. Because the only available 
experimental data for these states is for v — 0, we had to rely heavily on the ab initio 
calculations described above to supply the missing information on the vibrational structure. 
In particular, the first 4 vibrational levels were calculated from the ab initio potentials and 
fitted with 3 vibrational constants as reported at the end of Table 6. Note that the vibrational 
constants of Table 1 were generated in the same way but with 3 vibrational levels and 2 fitted 
vibrational constants (as listed in Table 1). The vibrational constants in Tables 1 and 6 are 
therefore slightly different. 

The electronic coupling in the X, A, and B states is intermediate Hund's case (a) - (b) 
over a wide range of J. The case (a) Hamiltonians for the 4 $ and 4 r states in Tables 3 and 4, 
respectively, were selected for least-squares fitting of the experimental data. In determining 
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the molecular constants for the X state, the reported rotational lines in Andersson et al. 
(2003) were reduced to lower state combination differences, yielding the molecular constants 
listed in Table 5. ^From these fitted constants, the X state term energies were calculated 
and then used in transforming the observed rotational lines in the A — X and 5-Ifl-O 
bands into term energies for the A and B states. Molecular constants derived from fits of 
the A and B term energies are also listed in Table 5. 

Except for the X state, the centrifugal constants D and H were not statistically 
determined. Instead, these constants were held fixed in the fits of the A and B states to the 
values determined by the Dunham approximations, 

D v ~ £ 2 B V and H v ~ ^uu e - fa e /3 (1) 

where £ = 2B v /u e and estimated values for u e and a e were supplied by our ab initio calcu- 
lations (Table 6). The presence of an anomalously large e/f parity splitting in the A state, 
undoubtedly due to perturbations, required separate fits of the e- and /-parity term energies, 
yielding the two sets of molecular constants in Table 5. Finally, the standard deviations from 
the fits, 0.031, 0.28, 0.16, and 1.18 cm" 1 for the X, A(e), A(f), and B states, are comparable 
in quality to 0.026 and 0.42 cm" 1 for X and A (Andersson et al. 2003) and 0.020 and 0.573 
cm -1 for X and B (Launila & Lingren 1996) with roughly the same number of data points 
and adjustable parameters. 

The molecular constants for v > were generated in a straightforward manner. With 
the ab initio estimates, cu e , u e x e , u> e y e , a e , j e , and e e , for each state (Table 6), the rotational 
constants for v > were computed from 

B v = B e -a e (v + 1/2) + le (v + l/2) 2 + e e (v + l/2) 3 (2) 

where the B e 's were determined from -Bo's, the centrifugal constants D v and H v from eq. (1), 
and the spin-component energies T V (Q) by starting from the T (Q) values and recursively 
adding vibrational separations of 

AG(v) = u e -2v uj e x e + (3v 2 + 1/4) uj e y e (3) 

for v > in which the dependence of A, A, and 7 on v was assumed to be negligible. The 
list of molecular constants in Table 6 were then used with the 4( J> and 4 r Hamiltonians in 
computing the term energies from v — to v — 5 for e/f rotational levels up to J = 50.5 
for the X, A(e), A(f), and B states. Note that all of the vibrational constants have been 
obtained using ab initio calculations so the location of bands with Av 7^ (computed using 
equation (3)) is relatively uncertain. 
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In our calculations the Einstein A-value is partitioned as 

A = A v i v ii HLF/(2J' + 1) 

where A v i v ii is the Einstein A for the vibronic transition v' — v", calculated using the ab 
initio electronic transition dipole moment function, and HLF is the Honl-London factor. 
For reference, the Honl-London factors (labeled as S(J)) in Tables 7-10 for 4 $ — 4 $ and 
4 r — 4 $ transitions involving pure Hund's case (a) and (b) coupling were derived using the 
method discussed in Dulick et al. (2003). Because the coupling in the X, A, and B states 
is intermediate Hund's case (a) - (b), the actual HLF's were computed by starting from 
U[TU U , where and U„ are the eigenvectors obtained from the diagonalizations of the 
lower and upper state Hamiltonians and T is the electric-dipole transition matrix (cf., eq. (1) 
in Dulick et al. (2003)). Squaring the matrix elements of U]TU U yield numbers proportional 
to the HLF's for intermediate coupling. In the absence of perturbations from neighboring 
states, as J becomes large, the intermediate case HLF's asymptotically approach the pure 
Hund's case (b) values. 



5. CHEMICAL ABUNDANCES OF TITANIUM HYDRIDE 

In order to estimate the contribution TiH makes to M dwarf spectra, its abundance must 
be calculated. This is accomplished by minimizing the total free energy of the ensemble of 
species using the same methods as Sharp & Huebner (1990) and Burrows & Sharp (1999) 
(see Appendix A) . To calculate the abundance of a particular species, its Gibbs free energy 
of formation from its constituent elements must be known. Using the new spectroscopic 
constants derived in §2 and §3, we have calculated for TiH this free energy and its associated 
partition function. In the case of TiH, its dissociation into its constituent elements can be 
written as the equilibrium 

TiH ^ Ti + H . 

The Gibbs energy of formation is then calculated from eq. (Al) in Appendix A. Eq. (Al) is 
the general formula for obtaining the energy of formation of a gas-phase atomic or molecular 
species, be it charged or neutral. 

Since the publication of Burrows & Sharp (1999), we have upgraded our chemical 
database significantly with the inclusion of several additional titanium-bearing compounds, 
and have updated most of the thermochemical data for titanium-bearing species (Barin 
1995). In addition to TiH, the new species are the solid condensates TiS, Mg 2 Ti0 4 , Ca 3 Ti 2 7 , 
and Ca 4 Ti 3 Oio, and the revised data are for the gas-phase species TiO, Ti0 2 , and TiS, the 
liquid condensates TiN, Ti0 2 , Ti 2 3 , Ti 3 5 , Ti 4 7 , MgTi0 3 , Mg 2 Ti0 4 , and the solid con- 
densates TiN, Ti0 2 , Ti 2 3 , Ti 3 5 , Ti 4 7 , MgTi0 3 , CaTi0 3 , and MgTi 2 5 . The revised 
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data for Ti0 2 , in particular, have resulted in refined abundance estimates of TiO. We have 
obtained much better least-square fits to the free energy, and in nearly all cases the fitted 
polynomials agree to better than 10 cal mole -1 compared with the tabulated values in Barin 
(1995) over the fitted temperature range, which for condensates corresponds to their stability 
field. For computational efficiency we sometimes extrapolate the polynomial fits of eq. (A2) 
for the condensates beyond the temperature range in which they are stable, but we ensure 
that the extrapolation is in a direction that further decreases stability. 

Since ionization plays a significant role at higher temperatures, a number of ionized 
species are included in the gas phase, including electrons as an additional "element" with a 
negative stoichiometric coefficient. The ions considered in our equilibrium code are e~, H + , 
H", Li+, Na+, K+ Cs+, Mg+, A1+, Si+, Ti+, Fe+, H+ Hj, OH+, CO+, NO+, N+, Nj, and 
H 3 + . Ions were not considered in the earlier work of Burrows and Sharp (1999). 

The partition function of TiH is obtained from the estimated spectroscopic constants 
in Tables 11 and 12 (data taken from Table 1 and Anglada et al. 1990, or simply guessed). 
Table 11 lists 15 electronic states, with the quartet X 4 $ ground state being resolved into each 
of the four separate spin substates. The partition function is calculated using the method 
of Dulick et al. (2003). The contribution to the total partition function of each electronic 
state is determined using asymptotic approximations from the vibrational and rotational 
constants in Table 12, then these contributions are summed according to eq. (7) in Dulick et 
al. (2003), with the Boltzmann factor for the electronic energy being applied to each state, 
where the electronic energy of the lowest spin substate of the ground X electronic state is 
zero by definition. Since the separate electronic energies for the spin substates of the ground 
electronic state are available, the X state is treated as four separate states in the summation, 
so that effectively the sum is over 18 states. The partition function of titanium is obtained 
by summing over the term values below 20,000 cm -1 , with degeneracies as listed in Moore 
(1949). The partition function of hydrogen is set to a value of 2 for all the temperatures 
considered. Table 13 lists our resulting TiH partition function in steps of 200 K from 1200 K 
to 4800 K. Using the partition functions of TiH, Ti, and H in eq. (Al), AG(T) is calculated 
at 100 K intervals over the temperature range required to compute the abundance of TiH, 
fitted to a polynomial, and then incorporated into our database. 

The thermochemistry indicates that the most important species of titanium are TiO, 
Ti0 2 , TiH, Ti, and Ti + . Ti0 2 is always subdominant and Ti and Ti + are important only 
at higher temperatures than obtain in M, L, or T dwarf atmospheres. At solar metallicity, 
just before the first titanium-bearing condensate appears, Ti0 2 replaces Ti as the second 
most abundant titanium-bearing species, reaching an abundance of ~9%. Figure 3 depicts 
the dependence of the TiO and TiH abundances with temperature from 1500 K to 5000 K 
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for a range of pressures (P) from 1CT 2 atmospheres to 10 2 atmospheres. This figure also 
shows the metallicity dependence of the TiO and TiH abundances. For solar metallicity, 
TiO dominates at low pressures below ~4000 K and at high pressures below ~4500 K, 
but the TiH abundances are not small. For solar metallicity, the TiH/TiO ratio at 10~ 2 
atmospheres and 2500 K is 10~ 4 and at 10 2 atmospheres and 2500 K it is ~10 -2 (one percent). 
The predominance of TiO before condensation, then its disappearance with condensation, 
matches the behavior seen in the M to L dwarf transition. 

However, the abundance of TiO is a more strongly decreasing function of temperature 
than is that of TiH. For a given pressure, there is a temperature above which the TiH/TiO 
ratio actually goes above one. As Fig. 3 indicates, this transition temperature decreases 
with decreasing metallicity. For 0.01 x solar metallicity and a pressure of 1 atmosphere, the 
TiH/TiO ratio hits unity at only ~3500 K and at 10 2 atmospheres it does so near 2500 
K. As a result, we expect that TiH will assume a more important role in low- metallicity 
subdwarfs than in "normal" dwarfs, with metallicities near solar. This fact is in keeping 
with the putative measurement of the TiH band near 0.94 //m (§6) in the subdwarf 2MASS 
J0532. 

As expected, the abundance of gas-phase species containing titanium drops rapidly when 
the first titanium-bearing condensates appear, which occurs for solar metallicity at tempera- 
tures between 2200 K and 1700 K for pressures of 10 2 and 10~ 2 atmospheres, respectively (see 
Fig. 3). At solar metallicity, when titanium-bearing condensates form, the TiH abundance 
is already falling with decreasing temperature; at condensation its abundance is at least an 
order of magnitude down from its peak abundance. As Fig. 3 indicates, the TiH abundance 
is an increasing function of pressure, but assumes a peak, sometimes broad, between the for- 
mation of condensates at lower temperatures and the formation of atomic Ti/Ti + at higher 
temperatures. 

The only significant previous study that included TiH is that of Hauschildt et al. (1997). 
These authors present numerical data in the form of fractional abundances in the gas phase 
for a number of titanium-bearing species, including TiO and TiH, as a function of optical 
depth. Unfortunately, this makes comparison with our results difficult. Nevertheless, we can 
still see that for their two solar-metallicity models the maximum ratio of the abundance of 
TiH to TiO is smaller than we obtain, sometimes by large factors. For instance, for their 
T e ff=2700 K model, at an optical depth of unity, their TiH abundance is over four orders 
of magnitude lower than that of TiO, and at an optical depth of 100, it is still down by 
three orders of magnitude. For their T e fj=4000 K model at unit optical depth, their TiH 
abundance is below that of TiO by three orders of magnitude. These very low abundances 
are in contrast to what we obtain at comparable T/P points (cf. Fig. 3). 
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On the other hand, we find at a temperature of 2700 K, a pressure of 100 atmospheres, 
and solar metallicity that the abundance of TiH is as high as ~1.5% of that of TiO. At 1 
atmosphere, this ratio has dropped by about a factor of 10, but is still about a factor of 
10 larger than the corresponding ratio in Hauschildt et al. At a temperature of 4000 K, a 
pressure of 100 atmospheres, and solar metallicity, we obtain a TiH/TiO abundance ratio as 
large as ~24%. Furthermore, as previously mentioned, we find that the TiH/TiO ratio can 
go above one and at low metallicities is above one for a wider temperature range. Even at 
temperatures where this ratio is significantly below one, we conclude that the TiH abundance 
derived with the new thermochemical data is much higher than previously estimated. 



6. TITANIUM HYDRIDE ABSORPTION OPACITIES 

Using the results of §2, §3, and §4, we calculate the opacity of the two electronic systems 
A 4 <3> - X 4 $ and B A T - X 4 $ of TiH. These have band systems centered near 0.94 /xmin 
the near IR and 0.53 //min the visible, respectively. For the A — X system, 4556 vibration- 
rotation transitions are considered for each of the 20 vibrational bands v' — v" with v'=0 to 
4 and v"=0 to 3, giving 91120 lines in total. For the B — X system, 4498 vibration-rotation 
transitions are incorporated for each of the 24 vibrational bands v'=0 to 3 and v"=0 to 5, 
giving 107,952 lines in total. 

We use the same methods to calculate the line strengths and broadening as in §5 of 
Dulick et al. (2003), using the partition function of TiH from §5. Since broadening param- 
eters for TiH are not available, those of FeH are used. In the absence of any other data, 
since TiH shares many characteristics with FeH, this is probably a reasonable assumption. 
A truncated Lorentzian profile is used, where the absorption in the wings of each line is 
truncated at min(25xP tot ,100) cm -1 on either side of the line center, P tot is the total gas 
pressure in atmospheres, and the detuning at which the absorption is set to zero does not 
exceed 100 cm -1 . To conserve the total line strength, the absorption in the part of the profile 
calculated is increased to account for the truncated wings. The reason for truncating the 
profiles is that we have no satisfactory theory to deal with the far wings of molecular lines, 
but we know that for atomic lines there is a rapid fall-off in the far wings, and, hence surmise 
that a rapid fall-off is probably realistic. However, our choice of determining its location is 
rather arbitrary. 

As with FeH, different isotopic versions are included in the calculations. The main 
isotope of Ti is 48 Ti, which makes up 73.8% of terrestrial titanium. The other isotopes of 
titanium are 46 Ti, 47 Ti, 49 Ti, and 50 Ti with isotopic abundances of 8%, 7.3%, 5.5%, and 
5.4%, respectively. Using the same methods as Dulick et al. (2003) based on Herzberg 
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(1950), the shift in the energy levels, and, hence, the transition frequencies for the minor 
isotopic versions of TiH are calculated. As in the case of Dulick et al., the change in the 
line strength with isotope depends only on the isotopic fraction of Ti in that particular 
isotopic form, and effects due to changes in the partition function, Boltzmann factor, and 
intrinsic line strength are small, and so are neglected. Figure 4 shows the absorption opacity 
in cm 2 per TiH molecule (plotted in red) between a wavelength of 0.3 /iinand 2.0 /xm, 
at temperature/pressure points of 3000 K/10 atmospheres and 2000 K/l atmosphere. For 
comparison, the H 2 (blue) opacities per molecule at 3000K and 10 atmospheres are also 
shown. The approximate location of the photometric bands B, V, and R in the visible, and 
J, J, and H in the infrared are also indicated. 

The absorption spectrum of TiH covering the region in the near infrared from ~0.7 
jum (the limit of the far visible red) to ~1.95 fxm , is due to the A — X electronic system, and 
the absorption between ~0.43 /zmand ~0.7 /mi in the visible is due to the B — X electronic 
system. The two electronic systems overlap at ~0.7 /im. The strongest absorption in the 
A — X system corresponds to the Av=0 bands (0-0, 1-1, 2-2, etc.), the next strongest peaks 
at a shorter wavelength and corresponds to the Av = 1 bands (1-0, 2-1, 3-2, etc.), and 
the next strongest peaks at a longer wavelength and corresponds to the Av = — 1 bands 
(0-1, 1-2, 2-3, etc.). There is a similar structure for the B — X system. Many of the peak 
absorption features of TiH occur at about the same wavelengths as the peaks in absorption 
of TiO, which is more abundant. This is particularly obvious for the A — X system of TiH. 
Moreover, although the cross section per molecule of H 2 is much less than that of TiH over 
most of the wavelength region, H 2 is far more abundant, and the five peaks corresponding 
to the Av =-2, -1, 0, 1 and 2 sequence of bands for this A — X system of TiH match up 
approximately with peaks in the H 2 absorption spectrum. 

Therefore, because of the particular wavelength positions of the TiO and H 2 absorption 
bands and the high abundances of the TiO and H 2 molecules at solar metallicity, identifying 
TiH in a solar-metallicity M or L dwarf spectrum, while not impossible, is likely to be difficult. 
However, the TiH features are not completely obscured, in particular at ~0.52 /mi, ~0.94 
/iiii, and weakly in the H band near 1.6 /im. Notice that the band heads of the A — X 
transition are mainly R-heads because B' < B" but for the Av = —3 sequence in the H 
band the shading changes and P-heads form. The peak in the TiH absorption in the H 
band corresponds with a deep trough of the TiO absorption, and is close to a minimum in 
the H 2 absorption. Therefore, in particular in subdwarfs with sub-solar metallicities, for 
which the abundances of metal hydrides come into their own vis a vis metal oxides (Fig. 
3), spectral features of TiH in the H band and at 0.94 /mi should emerge. The larger TiH 
abundances we derive in §5 also suggest this. The recent detection in the subdwarf L dwarf 
2MASS J0532 by Burgasser, Kirkpatrick & Lepine (2004) of the TiH feature at 0.94 /mi is 
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in keeping with these expectations. 

7. CONCLUSIONS 

For the TiH molecule, we have combined ab initio calculations with spectroscopic mea- 
surements to derive new thermochemical data, new spectral line lists and oscillator strengths, 
its abundance (for a given composition, temperature, and pressure), and its per-molecule ab- 
sorption opacities. We find that with the new partition functions, the abundance of TiH 
in M and L dwarf atmospheres is much higher than previously thought, that the TiH/TiO 
ratio increases strongly with decreasing metallicity, and that at high enough temperatures 
the TiH/TiO ratio can exceed unity. Furthermore, we conclude that, particularly for sub- 
dwarf L and M dwarfs, spectral features of TiH near 0.94 /zmand in the H band may not 
be weak. The recent putative detection of the 0.94 //m feature in the L subdwarf 2MASS 
J0532 is encouraging in this regard and suggests that the detection of TiH in other dwarfs 
may shed light on their atmospheric titanium chemistry. 
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Support was also provided by the NASA Laboratory Physics Program and the Natural Sci- 
ences and Engineering Research Council of Canada. The authors would like to thank R. 
Freedman for providing guidance with line broadening parameters. Furthermore, AB ac- 
knowledges support through the Cooperative Agreement #NNA04CC07A between the Uni- 
versity of Arizona/NOAO LAPLACE node and NASA's Astrobiology Institute. The line 
lists themselves are available in electronic form at http://bernath.uwaterloo.ca/TiH . 

A. FREE ENERGY OF A GAS-PHASE SPECIES 

The Gibbs free energy of formation of any gas-phase species, relative to its constituent 
element or elements in their monatomic gaseous state, can be calculated if its spectroscopic 
data, and those of its constituent elements, are known. We start with the equilibrium 

3 

Y« + qe- - J>^, 

i=i 

where Y represents an atomic or molecular species, and has the electric charge q, which 
is zero for a neutral species, Xi represents the i-th element in Y, rii is the corresponding 
stoichiometric coefficient, j is the number of different elements in Y, and e~ is the electron. 
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Based on a generalization of the Saha equation and the law of mass action, the following 
equation can be derived to obtain the Gibbs energy of formation of Y q : 



AG(Y q ) = -RT 



3 h r TP K 

In {Q[Y q \) "5> ln (Q[Xi\) + -TTfT + 5 (1 + Q - N) InT 



kT 



+ ^|ln riirriij - gm e | + q ^ln2 + ^ lnm e 

f 5 3 r 27r 1 1 

+(l + « + JV){ 5 lnt+ 5 ln —j -lnA } 



(Al) 



where Q[Y q ] and are the partition functions corresponding to Y q and the atoms X i: 

respectively, and E is the total energy (in wavenumbers) to fully dissociate a neutral molecule 
into its constituent atoms, if the species is a molecule. For a positive ion, the ionization 
potential (or, for a multiply charged positive ion, the sum of the ionization potentials) must 
be included in E. For a negative ion, the electron affinity must be subtracted from E. In 
addition, the total number of atoms in the molecule is N, and is given by ^21=1 n ii l - e -i the 
sum of the number of atoms of each element over the number of elements. The mass of the 
electron and the mass of the i-th element are m e and rrii in amu, respectively, R, Na, and A 
are the gas constant, Avogadro's number, and the standard atmosphere, respectively, and 
the other symbols have their usual meanings. 

To be rigorously correct, for ions the molecular weight is corrected for the gain/loss 
in mass due to the addition/removal of electrons. For neutral monatomic species all terms 
on the right hand side of eq. (Al) cancel yielding zero, and the energy of formation for 
free electrons is zero by definition. For convenience, all the temperature-independent terms 
(after multiplying by —RT) have been collected together, and depend only on the mass, the 
number of atoms, and charge of the species. 

In order to incorporate the free energies into our existing database in a uniform way, 
eq. (Al) is calculated at 100 K intervals over the temperature range we are likely to use, 
then a polynomial fit of the form 

AG(T) =a/T + b + cT + dT 2 + eT 3 (A2) 

is made. The lowest number of coefficients that make the best fit over the calculated range 
is adopted. The results for TiH using this procedure are a=-3.046x 10 5 , 6=-4.915x 10 4 , 
c=2.325xl0 1 , d=-1.568xl0~ 4 , and e=4.539xl0~ 8 for G(T) in calories per mole. 

In the case of diatomic neutral molecules, such as TiH, eq. (Al) simplifies by putting 
q = 0, N = 2, j = 2 and n\ = n 2 = 1, and the appropriate partition functions for the 



-16- 



molecule and the atoms are used. In addition, the energy E to fully dissociate the molecule 
is just the dissociation potential from the ground vibrational level of the ground electronic 
state, which for TiH is 2.08 eV, or 16776 cm" 1 . 
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Table 1: Summary of IC-MRCI+Q Spectroscopic Constants. 



State r (A) cj e (cm x ) uj e x e (cm x ) Dipole(debye) a T 00 (cm 1 ) 

valence correlation 









A 4 $ 


— X 4 $ calculation 








i 


.818 


1533.4 


21.94 


2.30 





A 4 $ 


i 


.896 


1375.9 
B 4 V 


23.17 
— X 4 $ calculation 


2.81 


10745 


x 4 $ 


i 


.821 


1527.0 


21.85 


2.13 





4 n 


i 


.838 


1489.3 


22.31 


1.91 


1694 


4 A 


i 


.905 


1402.4 


21.82 


1.17 


4082 


B 4 T 


i 


.804 


1577.8 


27.52 


3.19 


17735 








core+valence correlation 












A 4 $ 


— X 4 $ calculation 






X 4 $ 


i 


.788 


1548.9 


20.23 


2.46 





A 4 $ 


i 


.888 


1342.6 
B A Y 


21.26 
— X 4 $ calculation 


3.25 


11237 


X 4 $ 


i 


.794 


1529.2 


20.04 


1.52 6 





4 n 


i 


.813 


1487.4 


20.36 


1.40 


1732 


4 A 


i 


.885 


1398.7 


20.66 


1.24 


4231 


B 4 T 


i 


.764 


1592.5 


27.07 
Experiment 


2.59 


18874 


X 4 $ 


i 


.779 c 


1385.3 d 




2.455 e 





A 4 $ 


i 


.867 c 








10595' 


B 4 r 


i 


.7248' 






2.998 e 


18878' 



a Computed as an expectation value at the IC-MRCI level. 

b Using a finite field approach, the IC-MRCI(+Q) level yields 2.06(2. 17)debye. 

c Andersson et al. (2003). Note Launila & Lindgren (1996) give 1.7847 A for r of the X 

state. 

d Chertihin & Andrews (1994). 

e Steimle et al. (1991). 

f Launila & Lindgren (1996). 
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Table 2. Einstein A„/„« Values for the Bands of the B A T — X 4 $ and 

- X 4 $ Transitions. 



/ 

V 


If 

V 


A (s ) [BY — a $j 


A (s J [A <P — X <P J 


o 


o 


4003 xlO 8 


2336 xlO 7 

V 7 . '.J '.J V 7 /\ X V7 


o 


1 


6103xl0 6 


7567xl0 6 

V7 • 1 t_/ V7 1 s\ XV/ 


o 


2 


1279xl0 5 


1776 xlO 6 

Ui 1 1 ( V7 / \ XV/ 


n 


3 


n 9079x1 n 3 


n sqi Qx i n 5 


o 


4 


0.7146 




o 


5 


0.2469 




1 


o 


3239 xlO 6 


2577xl0 6 


1 


1 


3831 xlO 8 


1329xl0 7 

V 7 > X -.7 » 7 /\ X V 7 


1 


2 


1105xl0 7 


1071xl0 7 

\J • X \J 1 X /N X \J 


1 


3 


3751 xlO 5 


4074 xlO 6 

V7 . j: V7 1 i /\ X \J 


1 


4 


9085xl0 3 




1 


5 


6571X10 1 




2 





1205xl0 4 


8251 xlO 4 

V 7 > v7 ".7 X / \ X V 7 


2 


1 


5105 xlO 6 


4106xl0 6 


2 


2 


3664 xlO 8 


6651 xlO 6 


2 


3 


1463 xlO 7 

V7 • X \J '.J /\ X V7 


1090xl0 7 

\j • x vy y\ x w 


2 


4 


0.7186xl0 5 




2 


5 


0.2383xl0 4 




3 





0.2718xl0 3 


0.9007xl0 2 


3 


1 


0.8623xl0 4 


0.2273xl0 5 


3 


2 


0.5489xl0 6 


0.4793xl0 6 


3 


3 


0.3499xl0 8 


0.4793xl0 6 


3 


4 


0.1664xl0 7 




3 


5 


0.1119xl0 6 




4 







0.2572 


4 


1 




0.3647xl0 3 


4 


2 




0.4142xl0 5 


4 


3 




0.4816xl0 6 
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Table 3. 4 <£> Hund's Case (a) Hamiltonian Matrix 



x = J +1/2, t = ^/3(x 2 - 4), u = ^(a; 2 - 9), w = ^3(a; 2 - 16) 
T(9/2) = T + 9A/2 + 2A + 3 7 
T(7/2) = T + 3A/2 - 2A - 2 7 
T(5/2) = T - 3A/2 - 2A - 5 7 
T(3/2) = T - 9,4/2 + 2A - 67 



< 9/2 


H 


9/2 > 


= T(9/2) + (x 2 - 10) B - (x 4 


- I7x 2 + 52)D + (x 6 - 21x 4 + 96a; 2 - 40)H 


< 9/2 


H 


7/2 > 


= < 7/2|#|9/2 > = —v [B — 


7/2 - 2{x 2 - h)D + (3a; 4 - 23a; 2 + 16) H] 


< 9/2 


H 


5/2 > 


= <5/2|ff|9/2>= -uv[D 


- (3a; 2 - A)H] 


< 9/2 


H 


3/2 > 


= < 3/2|#|9/2 > = -tuvH 




< 7/2 


H 


7/2 > 


= T(7/2) + x 2 B-(x 4 + 7x 2 


- 84)£) + (a; 6 + 21a; 4 - 258a; 2 + 2M)H 


< 7/2 


H 


5/2 > 


= < 5/2|#|7/2 >= -u[B- 


7/2 - 2(a; 2 + 3)D + (3.x 4 + 28a; 2 - 60) H] 


< 7/2 


H 


3/2 > 


= < 3/2|#|7/2 > = -tu (D 


- (3a; 2 + U)H) 


< 5/2 


H 


5/2 > 


= T(5/2) + (x 2 + 6)B — (a; 4 - 


f 19a; 2 - 12)D + (.a; 6 + 39a; 4 + 72a; 2 - Ab6)H 


< 5/2 


H 


3/2 > 


= < 3/2|#|5/2 > = —t [B — 


7/2 - 2(a; 2 + 7)D + (3a; 4 + 49a; 2 + 100) H] 


< 3/2 


H 


3/2 > 


= T(3/2) + (x 2 + 8)B-(x 4 - 


f 19a; 2 + 52)L> + (a; 6 + 33a; 4 + 222a; 2 + 248)i? 
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Table 4. 4 r Hund's Case (a) Hamiltonian Matrix 



< ll/2|#|ll/2> 



< 11/2 

< 11/2 

< 11/2 

< 9/2 

< 9/2 

< 9/2 

< 7/2 

< 7/2 

< 5/2 



\H 
\H 
\H 
\H 
\H 
\H 
\H 
\H 
\H 



|9/2 > 
|7/2 > 
|5/2> 
|9/2> 
|7/2 > 
|5/2> 
|7/2> 
|5/2> 
|5/2> 



J +1/2, t = ^/3(x 2 - 25), u = ^(x 2 - 16), v = ^Ji{x 2 - 9) 

T(ll/2) = T + 6A + 2A + 9 7 /2 

T(9/2) = T + 2A-2A-3 7 /2 

T(7/2) = T - 2A - 2A - II7/2 

T(5/2) = T - 6A + 2A - 15 7 /2 

= T(ll/2) + (x 2 - 13)B - (a; 4 - 23a; 2 + 94)L> + (a; 6 - 30a; 4 + 201x 2 - 172)i? 

= < 9/2|ff|ll/2 > = -i [S - 7/2 - 2(x 2 - 7)D + (3a; 4 - 35a; 2 + 44) if] 

= < 7/2|if|ll/2 > = -tu [D - (3a; 2 - 7)if] 

= < 5/2|if|ll/2 > = -tuvH 

= T(9/2) + (a; 2 - 1)B - (a; 4 + 5a; 2 - 138)£> + (a; 6 + 18a; 4 - 439a; 2 + 804)if 

= < 7/2|JT|9/2 > = -u [B - 7/2 - 2(x 2 + 3)D + (3a; 4 + 28a; 2 - 123)if] 

= < 5/2|if |9/2 > = -uv [D - (3a; 2 + 17)H] 

= T(7/2) + (x 2 + 7)B - (x 4 + 21a; 2 - 42)L> + (a; 6 + 42a; 4 + a; 2 - 1164)if 

= < 5/2|if |7/2 > = —v [B — 7/2 - 2(x 2 + 9)D + (3a; 4 + 61a; 2 + 156) H] 

= T(5/2) + {x 2 + U)B - (a; 4 + 25a; 2 + 94)L> + (a; 6 + 42a; 4 + 369a; 2 + 548) if 
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Table 5. Derived Molecular Constants (cm 1 ) From Observed Data 



State 


To 


B 


D x 10 4 


H x 10 9 


A 


7o 


A 4 $(e) a 
B 4 T 


85.1362 
10688.3028(848) 
10687.880(156) 
18799.595(238) 


5.3668672(614) 
4.876910(846) 
4.86282(109) 
5.70361(115) 


2.61187(303) 
2.58010 
2.55780 
2.94312 


9.108(405) 
-8.025 
-8.112 
-7.459 


32.96454(797) 
30.4589(246) 
29.0992(410) 
39.1548(442) 


0.193432(889) 
-0.10251(733) 
-0.22281(435) 
-1.0047(147) 



a Spin-spin constants, A(e) = -1.599(109) and A(/) = -5.586(118). 
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Table 6. Generated Molecular Constants (cm -1 ) for the X 4 $, A 4 $, and B 4 r Vibrational 

Levels 



V 




B v 


D l 


, x 10 4 


H v x 10 9 


T„(A-3/2) T„(A-l/2) 


T„(A+l/2) 


T„(A + 3/2) 














j\. otate 











5, 


,3668672 


2 


.61187 


9.108 


-64.3649 


OA TOOO 


134.1962 


234.0570 


1 


■5. 


,1007701 


2 


.21602 


-8.378 


1/1/1/1 1 oil 
1444.1222 


1 c /i o onno 


1642.6833 


1742.5441 


2 


4. 


.8611877 


1 


.91820 


-8.781 


2912.148/ 


oni 1 ooco 
oUll.zooo 


3110.7098 


3210.5707 


3 


4. 


,6513671 


1 


.68039 


-8.813 


/I QQQ /I /I CO 

4oo8.448y 


44o < .DoOl) 


4537.0100 


4636.8709 


4 


4. 


.4745555 


1 


.49595 


-8.652 


r7ni 7C71 
1 21. / O 1 1 


coon o /i /i o 


5920.3182 


6020.1790 


■5 


4, 


,3340001 


1 


.35936 


-8.424 


7060.8075 


7159.8946 


7259.3686 


7359.2294 














A 4 $(e) State b 











4. 


.8769102 


2 


.58010 


-8.025 


10548.6548 


10646.3250 


10737.3941 


10821.8623 


1 


4, 


,6169393 


2 


.18919 


-9.213 


11848.7691 


11946.4393 


12037.5085 


12121.9766 


2 


4. 


,3736732 


1 


.86106 


-9.643 


13106.3657 


13204.0358 


13295.1050 


13379.5732 


3 


4, 


,1484432 


1 


.58810 


-9.585 


14319.7104 


14417.3805 


14508.4497 


14592.9179 


4 


3, 


,9425804 


1 


.36321 


-9.237 


15487.0692 


15584.7394 


15675.8086 


15760.2767 


5 


3, 


,7574161 


1 


.18002 


-8.745 


16606. /U82 


1 d*7C\ A O TO A 

16704.3784 


16795.4476 


16879.9157 














/I 4^/ ^\ n +nf „b 

/l sP(/ J otate 











4. 


,8628162 


2 


.55780 


-8.112 


inr/17 AOQO 


1 C\d X. 1 77 

11)000. Ol I ( 


10743.1469 


10806.9864 


1 


4, 


,6028453 


2 


.16921 


-9.254 


1 1 0/17 01 QQ 
1 184 / .Zloo 


1 lyoO.Oozl 


12043.2612 


12107.1008 


2 


4. 


,3595792 


1 


.84313 


-9.651 


1 Q 1 H/1 QflOQ 

lolU4.oUyo 


1 QOI A OOQA 
Iozl4.zz80 


13300.8578 


13364.6973 


3 


4, 


.1343492 


1 


.57197 


-9.569 


14318.1545 


14427.5733 


14514.2025 


14578.0420 


4 


3 


,9284864 


1 


.34864 


-9.205 


15485.5134 


15594.9321 


15681.5613 


15745.4009 


5 


3, 


,7433221 


1 


.16679 


-8.701 


16605.1524 


16714.5711 


16801.2003 


16865.0399 














B 4 r State c 











5, 


,7036054 


2 


.94312 


-7.459 


18572.2017 


18726.8115 


18879.4119 


19030.0029 


1 


5. 


,4171527 


2 


.52009 


-8.938 


20110.5248 


20265.1346 


20417.7350 


20568.3260 


2 


5, 


,1537792 


2 


.17010 


-9.613 


21594.7116 


21749.3214 


21901.9218 


22052.5128 


3 


4, 


,9169618 


1 


.88449 


-9.772 


23020.5977 


23175.2075 


23327.8079 


23478.3989 


4 


4, 


,7101772 


1 


.65659 


-9.636 


24384.0185 


24538.6283 


24691.2287 


24841.8197 


5 


4. 


,5369023 


1 


.48041 


-9.363 


25680.8094 


25835.4192 


25988.0196 


26138.6106 



a w e = 1547.7347, u e x e = 19.2810, Lo e y e = -0.210954, B e = 5.5088441, a e = 0.2895, 7e = 0.010822, 
e e = 0.00054119, and 70-5 = 0.193432. 
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V = 1340.9704, u e x e = 19.9584, uj e y e = -0.289002, B e (e) = 5.012744, B e (f) = 4.99865, a e = 
0.2754, 7e = 0.007354, e e = 0.00022187, 70-5 (e) = -0.10251, and 70-5 (/) = -0.22281. 

c uj e = 1588.4683, u e x e = 23.9447, u e y e = -0.694092, B e = 5.8544, a e = 0.3062, 7e = 0.008932, 
e e = 0.00057947, 70-5 = -1.0047. 
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Table 7. 4 $ - 4 $ Honl-London Factors - Case (a) Coupling 



Branch 


fl" 


Parity 


0' 


Parity 




S(J) 


P(J) 


9/2 


e/f 


9/2 


e/f 


1 


(2J-9) (2J + 9) 


4 ' 


J 


P(J) 


7/2 


e/f 


7/2 


e/f 


1 

4 ' 


(2 J - 7) (2J + 7) 
J 


P(J) 


5/2 


e/f 


5/2 


e/f 


1 

4 ' 


(2 J - 5) (2J + 5) 
J 


P(J) 


3/2 


e/f 


3/2 


e/f 


1 

4 ' 


(2J-3) (2J + 3) 
J 


Q(J) 


9/2 


e/f 


9/2 


f/e 


81 
T 


2J+ 1 
J(J + 1) 


Q(J) 


7/2 


e/f 


7/2 


f/e 


49 

T 


2J+ 1 
J(J + 1) 


Q(J) 


5/2 


e/f 


5/2 


f/e 


25 

T 


2J+ 1 
J(J+1) 


Q(J) 


3/2 


e/f 


3/2 


f/e 


9 

4 ' 


2J+ 1 
J(J+1) 


R(J) 


9/2 


e/f 


9/2 


e/f 


1 

4 ' 


(2J-7) (2J+11) 
J+ 1 


R(J) 


7/2 


e/f 


7/2 


e/f 


1 


(2J-5) (2J + 9) 


4 ' 


J+l 


R(J) 


5/2 


e/f 


5/2 


e/f 


1 

4 ' 


(2J-3) (2J + 7) 
J+l 


R(J) 


3/2 


e/f 


3/2 


e/f 


1 (2J- 1) (2J + 5) 


4 


J+l 
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Table 8. 4 $ - 4 $ Honl-London Factors - Case (b) Coupling 



Branch 


N" 


Parity 


N' 


Parity 


P P{J) 


J -3/2 


+/- 


J -5/2 


+/- 


Q P{J) 


J -3/2 


+/" 


J - 3/2 


+/- 


R P{J) 


J - 3/2 


+/" 


J - 1/2 


+/" 


P P(J) 


J - 1/2 


+/- 


J - 3/2 


+/- 


Q P(J) 


J - 1/2 


+/- 


J - 1/2 


+/- 


R P(J) 


J - 1/2 


+/" 


J +1/2 


+/" 


P P{J) 


J +1/2 


+/- 


J - 1/2 


+/- 


Qp{J) 


J +1/2 


+/" 


J +1/2 


+/" 


P P{J) 


J + 3/2 


+/- 


J + 1/2 


+/- 


Q Q(J) 


J - 3/2 


+/" 


J - 3/2 


+/" 


R Q(J) 


J - 3/2 


+/- 


J - 1/2 


+/- 


P Q(J) 


J - 1/2 


+/- 


J - 3/2 


+/- 


Q Q{J) 


J - 1/2 


+/" 


J - 1/2 


+/" 


R Q(J) 


J - 1/2 


+/- 


J +1/2 


+/- 


P Q(J) 


J +1/2 


+/- 


J - 1/2 


+/- 


Q Q(J) 


J +1/2 


+/" 


J +1/2 


+/" 



1 

4 

108- 

3 
4 
1 

4 

576- 

3 
4 
1 

4 

108 
1 

4 ' 
36 • 
3 

4 ' 
3 

4 

36- 
1 

4 ' 
1 

4 

36- 



s_(J) 

(2J-9) (2J+1) (2J + 3) 
(2J-3) (J - 1) 
2J+ 1 
J(2J-3) (2J- l) 2 
(2J-7) (2J + 5) 
J 2 (2J- l) 2 (J - 1) 

(2 J - 7) (2 J - 3) (2 J + 1) (2 J + 5) ( J + 1) 
J 2 (2J- l) 2 
(J-1)(J + 1) 



J(2J- l) 2 (2J+ l) 2 
(2J-5) (2J + 7) 



J 2 (2J+1) 2 (J+l) 

(2 J - 5) (2 J - 1) (2 J + 3) (2 J + 7) ( J - 1) 



J 2 (2J+1) 2 



2J- 1 



J(2J+1) 2 (2J + 3) 
(2J-3) (2J- 1) (2J + 9) 

(2J + 3) (J + l) 
(2J+1) (J + l) 



J(2J-1) 2 
(2J-7) (2J + 



l)(2J + 5) 



J 2 (2J- l) 2 
(2J-7) (2J+l)(2J + 5) 
J 2 (2J- l) 2 
(2J 2 + J-4) 2 
J(2J- l) 2 (2J+ 1) (J+l) 
(2J- 5) (2J- 1) (2J + 3) (2 J 



7) 



J 2 (2J+1)(J+1) 2 
(2 J - 5) (2 J - 1) (2 J + 3) (2 J + 7) 
J 2 (2J+1)(J+1) 2 
(2J 2 + 3 J - 3) 2 
J(2J+1) (2J + 3) 2 (J+l) 
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Table 8 — Continued 



Branch N" Parity N' Parity 



S(J) 



R Q(J) 


J +1/2 


+/" 


J + 3/2 


+/ 


P Q(J) 


J + 3/2 


+/- 


J +1/2 


+/ 


Q Q(J) 


J + 3/2 


+/- 


J + 3/2 


+/ 


R R{J) 


J - 3/2 


+/- 


J - 1/2 


+/ 




J - 1/2 


+/- 


J - 1/2 


+/ 


R R(J) 


J- 1/2 


+/- 


J +1/2 


+/ 


P R{J) 


J + 1/2 


+/- 


J - 1/2 


+/ 


Qr(j) 


J +1/2 


+/- 


J +1/2 


+/ 


R R(J) 


J +1/2 


+/- 


J + 3/2 


+/ 


P R(J) 


J + 3/2 


+/- 


J +1/2 


+/ 


Qr(j) 


J + 3/2 


+/- 


J + 3/2 


+/ 


R R(J) 


J + 3/2 


+/- 


J + 5/2 


+/ 



3 

4 ' 
3 

4 ' 
36 
1 

4 ' 



(2J-3) (2J+1) (2J + 9) 

(2J + 3) 2 (^ + l) 2 
(2J-3) (2J+1) (2J + 9) 



(2J + 3) 2 (J ■ 
J(2J + 1) 



1)2 



(2J + 3) 2 (J + l) 
(2J-7) (2J + 3) (2J + 5) 



108- 



J(2J- 1) 
2J + 3 



1 

4 
3 
4 



(2J-1)(2J+1) 2 (J + 1) 
(2 J - 5) (2 J - 1) (2 J + 3) (2 J + 7) ( J + 2) 



(2 J 



(2 J +1) 2 (J +1)2 

5) (2J + 7) 



576- 



J(2J+1) 2 (J+1)2 

J (J + 2) 



1 

4 ' 
3 

4 ' 
108 



(2J+l)2(2J + 3) 2 (J+l) 
J (2 J - 3) (2 J + 1) (2 J + 5) (2 J - 



9) 



(2J + 3) 2 (J +1)2 
(2J-3) (2J + 9) 
(2J + 3) 2 (J + 1)2 (J + 2) 
2J+ 1 

' (J + l) (2J + 3) 2 (2J + 5) 
1 (2J-1) (2J+1)(2J+11) 
4' (2J + 5)(J + 2) 
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Table 9. 4 r — 4 $ Honl-London Factors - Case (a) Coupling 



Branch 





Parity 




Parity 


P(J) 


9/2 


e// 


11/2 


e/f 


P(J) 


7/2 


e// 


9/2 


e/f 


P(J) 


5/2 


e// 


7/2 


e/f 


P(J) 


3/2 


e// 


5/2 


e/f 


Q(J) 


9/2 


e// 


11/2 


f/e 


Q(J) 


7/2 


e/f 


9/2 


f/e 


Q(J) 


5/2 


e/f 


7/2 


f/e 


Q(J) 


3/2 


e/f 


5/2 


f/e 


R(J) 


9/2 


e/f 


11/2 


e/f 


R(J) 


7/2 


e/f 


9/2 


e/f 


R(J) 


5/2 


e/f 


7/2 


e/f 


R(J) 


3/2 


e/f 


5/2 


e/f 



S(J) 



1 


(2J 


- 11) (2J-9) 


8 




J 


1 


(2J 


-9) (2J-7) 


8 




J 


1 


(2J 


- 7) (2 J - 5) 


8 




J 


1 


(2J 


-5) (2J-3) 


8 




J 


1 


(2J 


- 9) (2J+ 1) (2J+ 11) 


8 




J(J+1) 


1 


(2J 


-7) (2J+1) (2J + 9) 


8 




J(J+1) 


1 


(2J 


-5) (2J+1) (2J + 7) 


8 




J(J+1) 


1 


(2J 


-3) (2J+1) (2J + 5) 


8 




J(J+1) 


1 


(2J + 11) (2J+13) 


8 




J+l 


1 


(2J 


+ 9) (2J+11) 


8 




J+l 


1 


(2J 


+ 7) (2J + 9) 


8 




J+l 


1 


(2J 


+ 5) (2J + 7) 


8 




J+l 
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Table 10. 4 r - 4 $ H6nl-London Factors - Case (b) Coupling 



Branch 


N" 


Parity 


iV' 


Parity 


P P{J) 


J -3/2 


+/- 


J - 5/2 


+/- 


Qp(J) 


J - 3/2 


+/" 


J - 3/2 


+/" 


R P{J) 


J -3/2 


+/- 


J - 1/2 


+/- 


p p(j) 


J - 1/2 


+/- 


J - 3/2 


+/" 


Qp(J) 


J - 1/2 


+/" 


J - 1/2 


+/- 


R P(J) 


J - 1/2 


+/" 


J + 1/2 


+/- 


P P(J) 


J + 1/2 


+/- 


J - 1/2 


+/" 


Qp(J) 


J +1/2 


+/- 


J +1/2 


+/- 


Pp(j) 


J + 3/2 


+/" 


J +1/2 


+/- 


Q Q{J) 


J - 3/2 


+/- 


J - 3/2 


+/" 


R Q(J) 


J - 3/2 


+/" 


J - 1/2 


+/" 


P Q(J) 


J - 1/2 


+/- 


J - 3/2 


+/" 


Q Q(J) 


J - 1/2 


+/" 


J - 1/2 


+/" 


R Q(J) 


J - 1/2 


+/- 


J +1/2 


+/- 


P Q(J) 


J + 1/2 


+/" 


J - 1/2 


+/" 


Q Q(J) 


J + 1/2 


+/" 


J + 1/2 


+/" 



1 (2J- 11) (2J-9) (2J+1) 



(2J-3) (J - 1) 
(2J-9) (2J+l)(2J + 5) 
J(2J-3) (2J- l) 2 
(2J+5) (2J + 7) 
' J 2 (2J-1) 2 (J-l) 
(2 J - 9) (2 J - 7) (2 J - 3) (2 J + 1) (J + 1) 

J 2 (2J- l) 2 
(2J-7) (2J + 7) (J-l) (J + l) 
J(2J-1) 2 (2J+1) 2 
(2J + 7) (2J + 9) 
' J 2 (2J+1) 2 (J+l) 
(2 J - 7) (2 J - 5) (2 J - 1) (2 J + 3) ( J - 1) 



(2J- 5) (2J- 



J 2 (2J+ l) 2 
1) (2J + 9) 



J(2J+ l) 2 (2J + 3) 
(2J-5) (2J-3) (2J- 1) 

(2J + 3) (J + l) 
(2 J - 9) (2 J + 5) (2 J + 1) ( J + 1) 

J(2J- l) 2 
(2J+ 1) (2J + 5) (2J+ 7) 

J 2 (2J-1) 2 
(2J-9) (2J-7) (2J+1) 

J 2 (2J- l) 2 
(2J-7) (2J 2 + J-4) 2 (2J + 7) 

J(2J- l) 2 (2J+ 1) (J+l) 
(2 J - 1) (2 J + 3) (2 J + 7) (2 J + 9) 

J 2 (2J+1)(J+1) 2 
(2 J - 7) (2 J - 5) (2 J - 1) (2 J + 3) 



(2J- 



J 2 (2J+1)(J+1) 2 
5) (2 J 2 + 3J- 3) 2 (2J + 9) 



J(2J+1)(2J + 3) 2 (J+1) 
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Table 10 — Continued 



Branch 


N" 


Parity 


TV' 


Parity 




5(J) 


R Q(J) 


J +1/2 


+/" 


J + 3/2 


+/" 


3 


(2J+ 1) (2J + 9) (2J+ 11) 


8 


(2J + 3) 2 (J+ l) 2 


P Q(J) 


J + 3/2 


+/" 


J +1/2 


+/" 


3 


(2J-5) (2J-3) (2J+ 1) 


8 


(2J + 3) 2 (J + l) 2 


Q Q{J) 


J + 3/2 


+/" 


J + 3/2 


+/" 


1 

2 


J (2 J - 3) (2 J + 1) (2 J + 11) 
(2J + 3) 2 (J+l) 


R R(J) 


J -3/2 


+/- 


J - 1/2 


+/- 


1 


(2J + 3) (2J + 5) (2J + 7) 


8 


J(2J- 1) 




J - 1/2 


+/" 


J - 1/2 


+/- 


3 
2 


(2J- 7) (2J + 3) (2J+7) 
(2J-1) (2J+1) 2 (J+1) 


R R(J) 


J- 1/2 


+/" 


J +1/2 


+/- 


1 


(2 J - 1) (2 J + 3) (2 J + 7) (2 J + 9) ( J + 2) 


8 


(2J+1) 2 (J+l) 2 


P R{J) 


J + 1/2 


+/- 


J - 1/2 


+/" 


3 


(2 J - 7) (2 J - 5) 


8 


J(2J+1) 2 (J + 1) 2 


Qr(j) 


J +1/2 


+/" 


J +1/2 


+/" 


8 • 


J(2J-5) (2J + 9) (J + 2) 
(2J+1) 2 (2J + 3) 2 (J+l) 


R R(J) 


J +1/2 


+/- 


J + 3/2 


+/- 


1 


J (2 J + 1) (2J + 5) (2J + 9) (2J + 11) 


8 


(2J + 3) 2 (J + l) 2 


P R{J) 


J + 3/2 


+/- 


J +1/2 


+/" 


3 


(2J-5) (2J-3) 


8 


(2J + 3) 2 (J + l) 2 (J + 2) 


Qr(j) 


J + 3/2 


+/" 


J + 3/2 


+/- 


3 


(2J-3) (2J+ 1) (2J+ 11) 


2 


(J + l) (2J + 3) 2 (2J + 5) 


R R{J) 


J + 3/2 


+/" 


J + 5/2 


+/" 


1 

8 


(2J+ 1) (2J+ 11) (2J+ 13) 
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Table 11. Estimated TiH Data for 
Thermo chemistry 



State Term value in cm 



X 4 $ 3/2 


0.0 


X 4 $5/2 


99.6 


x 4 $ 7/2 


199.2 


X 4 $ 9/2 


298.8 


1 4 S~ 


970 


i 4 n 


1732 


1 2 A 


1770 


1 4 A 


4231 


i 2 n 


3390 


i 2 $ 


3710 


i 2 s- 


5810 


2 4 £~ 


8630 


2 4 A 


9520 


2 4 n 


10810 


a 4 $ 


10595 


2 2 n 


6940 


2 2 A 


13870 


B 4 T 


18878 
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Table 12. Estimated Molecular Constants of TiH for 
Thermochemistry (in cm -1 ) 



State 




a e 






u e x e 




5.4975 


0.286 


2.56xl0- 4 


1548.9 


20.2 


1 4 S~ 


5.08 


0.27 


2.5xl0~ 4 


1484 


22.8 


i 4 n 


5.32 


0.29 


2.5xl0~ 4 


1487 


20.4 


1 2 A 


5.36 


0.29 


2.5xl0~ 4 


1519 


26.9 


1 4 A 


4.92 


0.27 


2.6xl0~ 4 


1399 


20.6 


i 2 n 


5.22 


0.28 


2.5xl0~ 4 


1450 


24.6 


i 2 $ 


5.00 


0.27 


2.5xl0~ 4 


1462 


26.6 


1 2 £- 


5.02 


0.27 


2.5xl0~ 4 


1446 


27.4 


2 4 S- 


5.13 


0.28 


2.5xl0~ 4 


1469 


27.4 


2 4 A 


4.81 


0.27 


2.5xl0~ 4 


1367 


23.2 


2 4 n 


4.86 


0.27 


2.5xl0~ 4 


1403 


29.6 


a 4 $ 


4.9915 


0.274 


2.6xl0~ 4 


1342.6 


21.3 


2 2 n 


5.25 


0.29 


2.5xl0~ 4 


1564 


24.7 


2 2 A 


5.03 


0.27 


2.5xl0~ 4 


1435 


22.1 


b 4 t 


5.8906 


0.303 


3.6xl0~ 4 


1592.5 


27.1 
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Table 13. Derived Partition Function of TiH 



lemperaxure 


rarxixion r unction 


1200 


1852 


1400 


2560 


1600 


3427 


1800 


4469 


2000 


5700 


2200 


7136 


2400 


8792 


2600 


10682 


2800 


12825 


3000 


15236 


3200 


17930 


3400 


20926 


3600 


24240 


3800 


27890 


4000 


31895 


4200 


36272 


4400 


41040 


4600 


46217 


4800 


51824 
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r(Ti-H) (A) 

Fig. 1. — The core+valence IC-MRCI+Q potential energy curves for the X, A, and B states. 
The zero of energy is Ti and H at infinite separation. The X state is plotted twice, one curve 
is from the A — X calculations while the second is from the B-X calculations. See the text 
for a full description of the two sets of calculations. 



-37- 




Fig. 2.- 



The core+valence IC-MRCI A — X and B — X transition dipole moments. 
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Fig. 3. — The log (base 10) of the number fractions of both TiH (solid) and TiO (dashed) 
as a function of temperature (T) from 1500 to 5000 Kelvin, for pressures from 10~ 2 at- 
mospheres to 10 2 atmospheres in equal logarithmic steps. Pressure increases for a given 
sequence from bottom to top. Furthermore, the condensation lines march rightward as the 
pressure increases. The abundances of both species are given for solar metallicity (blue), 
0.1 x solar metallicity (green), and 0.01 x solar metallicity (red). TiO generally dominates at 
lower temperatures, until titanium condensables, such as perovskite, form, after which the 
abundances of all titanium compounds decrease precipitously (left of diagram). The con- 
densation temperature increases with increasing pressure. However, at high temperatures, 
TiH comes into its own and its abundance can exceed that of TiO. For a given pressure, 
the temperature above which the TiH/TiO ratio exceeds unity is an increasing function of 
metallicity, being at 10 2 atmospheres and solar metallicity ~4500 K, but at 10 2 atmospheres 
and 0.01 x solar metallicity ~2500 K. As metallicity decreases, TiH becomes more and more 
important, though the abundances of both TiO and TiH decrease with metallicity. When 
the first titanium-bearing condensate appears, the abundance of TiH is already falling. Note 
that with increasing pressure the peak in TiH abundance shifts to higher temperatures more 
rapidly than the condensation temperature. 
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Fig. 4. — The log (base 10) of a, the absorption cross section per molecule in cm 2 , for TiH 
at 10 atmospheres and 3000 K (red) and 1 atmosphere and 2000 K (black). For comparison, 
the cross section per molecule for H2O at 10 atmospheres and 3000 K is shown in blue. The 
approximate positions of the photometric bands B, V, R, I, J, and H are also indicated. 
The broad absorption of TiH with several peaks between ~0.7 /imand ~1.95 /imis due to 
the A — X electronic band system, and the absorption between ~0.43 /imand ~0.7 /imis 
due to the B — X electronic band system. See text for a detailed discussion. 



